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■ Abstract 

(— | . Spatially homogeneous but totally anisotropic and non-flat Bianchi type II cosmological model has been 

studied in general relativity in the presence of two minimally interacting fluids; a perfect fluid as the matter fluid 
and a hypothetical anisotropic fluid as the dark energy fluid. The Einstein's field equations have been solved by 

■ applying two kinematical ansatze: we have assumed the variation law for the mean Hubble parameter that yields 
a constant value of deceleration parameter, and one of the components of the shear tensor has been considered 
proportional to the mean Hubble parameter. We have particularly dwelled on the accelerating models with 
non-divergent expansion anisotropy as the Universe evolves. Yielding anisotropic pressure, the fluid we consider 
in the context of dark energy, can produce results that can be produced in the presence of isotropic fluid in 
accordance with the ACDM cosmology. However, the derived model gives additional opportunities by being able 
to allow kinematics that cannot be produced in the presence of fluids that yield only isotropic pressure. We have 
obtained well behaving cases where the anisotropy of the expansion and the anisotropy of the fluid converge 
to finite values (include zero) in the late Universe. We have also showed that although the metric we consider 
is totally anisotropic, the anisotropy of the dark energy is constrained to be axially symmetric, as long as the 
overall energy momentum tensor possesses zero shear stress. 

00. 
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1 Introduction 

^—1 \ The discovery of the current acceleration of the Universe from the Type la Supernova (SNIa) observations [TJ [2] 
i— I ■ posed one of the most challenging problems in modern cosmology. This acceleration of the Universe is, today, a well 
established fact that is confirmed by various independent observational data, including SNIa, CMB radiation, etc. 
• <-h ■ However, this discovery can be maintained in the context of general relativity (GR) by introducing a mysterious 
kind of energy source called the dark energy (DE), that can generate repulsive gravity (see references [31 H] for 
extensive reviews on DE). It is well known that perfect fluids with a constant equation of state (EoS) parameter 
lower than — -| are, today, usually considered in the context of DE in cosmology, since they can give rise to an 
accelerating expansion. Beside such perfect fluids with constant EoS parameter, various scalar field models, which 
can conveniently be described by an isotropic and time dependent EoS parameter, have been studied with an 
intention of producing kinematics consistent with the observed acceleration of the Universe [3j 0]. However, this is 
not the only possible method of studying DE. It can sometimes be useful to use a reverse method, i.e., one may 
seek for the solutions of the Einstein field equations by introducing some kinematical ansatze that are consistent 
with the observed kinematics of the Universe, and may then obtain and investigate the dynamics of the fluid as a 
possible candidate of the DE. For instance, as an ansatz in the context of latter method, Berman's law [5j |6] for 
the Hubble parameter, that yields constant deceleration parameter (DP), has been widely considered for obtaining 
accelerating cosmological models explicitly in the framework of spatially homogenous spacetimes after the discovery 
of the current acceleration of the Universe (e.g., see Refs. [7HHH2] for recent studies on constant DP). 

The sources with anisotropic stresses have been rarely studied in cosmology until recent years. This might be 
because of the expectation that a dominance of an anisotropic stress would give rise to an anisotropic expansion 
contrary to observations which favor an Universe that expands isotropically. One may see references [101 E] for the 
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effects of some known anisotropic stresses (such as magnetic fields, collisionless relativistic particles, hydrodynamic 
shear viscosity etc.) on cosmological evolution. On the other hand, the recent developments in cosmology increased 
the interest in anisotropic stresses particularly in the context of DE. Studies focusing on DE, including those in the 
framework of spatially homogenous but anisotropic spacetimes, have generally assumed that DE yields an isotropic 
pressure and obeys a simple EoS in the form p — wp, where p is the energy density, p is the isotropic pressure and 
w is the EoS parameter, which is not necessarily constant. However, various discussions [TH [T31 Q3J [TSJ [THl E] that 
focused on whether the Wilkinson Microwave Anisotropy Probe (WMAP) data [TSJ HH [501 121] for more successful 
explanation need a Bianchi type morphology rather than Robertson- Walker (RW) type, have promoted the general 
interest not only in the Bianchi type cosmological models but also in the possibility of anisotropic nature of the 
DE. For instance, recently, Campanelli et al. |H] investigated the possibility that a DE component with anisotropic 
EoS could generate an ellipsoidal Universe with the correct characteristics to explain the low quadrupole in the 
CMB fluctuations and proposed that an ellipsoidal Universe allowed better matches with the large-scale CMB 
anisotropy data. This conclusion is also consistent with their more recent study [23j . where they analyzed the 
magnitude-redshift data of SNIa included in the Union catalogue by considering Bianchi type I spacetime in the 
presence of a DE fluid with anisotropic EoS. In this study they concluded that supernova data are compatible 
with a standard isotropic Universe, however a large level of anisotropy, both in the geometry of the Universe and 
in the EoS of DE, is allowed. Most recently, while this paper is in progress, Axelsson et al. [21] studied possible 
imperfect nature of DE and the large-angle anomalous features in the CMB by considering the 7-year WMAP 
temperature observations data. Their results also motivate the further investigation on the possible origin and 
constraints of imperfect source terms including imperfect DE fluid in cosmology. Following such studies, several 
cosmological models ESJ ESI ETJ [28l [29j [30] that introduce anisotropic DE candidates to break the spatial isotropy 
of the Universe, which is generally believed to be obtained via inflation in the early Universe, have been proposed. 
On the other hand, even if it turns out that, in the future, the RW spacetimes are exactly consistent with the 
observational data, this does not necessarily rule out the possible anisotropic nature of the DE. The reason being 
that anisotropic energy sources do not necessarily promote anisotropy in the expansion, and sometimes can even 
support isotropization |31[ 132) . 

Maximal spatial symmetry of the RW metrics does not allow the theoretical study of cosmological models 
containing anisotropic fluids as long as fluids are comoving. We indeed do not expect the fluids to have bulk 
motion at cosmological scales, and therefore we must take the spatially anisotropic spacetimes into account to be 
able to study cosmological anisotropic energy-momentum sources such as magnetic field, cosmic strings, as well 
as hypothetical fluids considered in the context of DE. In the context of the latter method we mentioned above, 
one may allow the fluid to be general as much as the metric allows and then seek for the solutions by introducing 
kincmatical ansatze for the evolution of the Universe rather than constraining the fluid from the beginning, for 
instance, by assuming that the fluid is isotropic. To obtain explicit cosmological models in GR, we do not have a 
general metric for spatially homogenous but anisotropic spacetime. Instead, we have nine of the Bianchi type metrics 
plus Kanstowski-Sach metric that form the complete class of spatially homogenous but not necessarily isotropic 
four-dimensional spacetimes [33, 34 . Hence, choosing one of these ten metrics as the geometry of the Universe gives 
opportunity to generalize the perfect fluid representation of the energy- momentum tensor (EMT). However, the 
symmetry properties of the chosen metric would still constrain EMTs that are allowed. Therefore, it is important to 
study cosmological models within the framework of different types of spatially anisotropic spacetimes, for a better 
understanding of the possibly anisotropic nature of the DE and its effects on the evolution of the Universe. Exact 
general relativistic cosmological models in the presence of anisotropic fluid, in the context of DE, using kinematical 
ansatze have been studied in the framework of the following metrics: Bianchi I 29, 30, 35, 36, 37, 38 , Bianchi VIo 
[391 E01 SU S3 S3], Bianchi III (31j Hfl [45] , Bianchi V EH [47], Kantowski- Sachs [HI HH [50], Bianchi VI [32] and 
locally rotationally symmetric (LRS) Bianchi type II [ST]. Such models have also been studied in f(R) theory of 
gravity instead of GR within the framework of Bianchi VIo [52] and Bianchi III [53] metrics. 

In this paper, we present general relativistic cosmological models within the framework of spatially homogeneous 
but totally anisotropic and non-flat Bianchi type II spacetime in the presence of two fluids, which are the minimally 
interacting perfect fluid that can be described by an isotropic and constant EoS parameter, and the dynamically 
anisotropic hypothetical fluid as the DE fluid. We make use of two kinematical ansatze at the same time, to fully 
determine the Einstein's field equations. As the first kinematical ansatz, we consider the generalized Berman's law 
[5] [6] for the Hubble parameter, that yields constant mean DP and whose negative values generate accelerating 
volumetric expansion. A similar procedure was carried out by Singh and Kumar [54], within the framework of 
LRS Bianchi II spacetime, in the presence of a single isotropic fluid with a constant EoS parameter. In a recent 
paper, Naidu et al. [SSJ have studied an LRS Bianchi type II DE model in a scalar tensor theory with generalized 
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Berman's law and variable EoS parameter of DE. We would also like to note that although the constant DP is not 
capable of representing the whole or a long period of history of the Universe at once, it can still be considered as 
an approximation for representing a particular period of the history of the Universe, for instance, the vicinity of 
the present time of the Universe. As the second kinematical ansatz, we consider the assumption first introduced by 
Bali and Jain [56] , that assumes one of the component of the shear tensor to be proportional to the mean Hubble 
parameter. 



2 The metric, field equations, solution and some immediate results 

We assume that the GR is valid at cosmological scales and consider Einstein's field equations 



8ttG. 
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where G\ is the Einstein tensor, G is the Newton's gravitational constant, c is the speed of light in the vacuum, R\ 
is the Ricci tensor, R — R\ is the Ricci scalar, 8f is the unit four-tensor and T? is the overall EMT that describes 
the physical ingredients of the Universe. 

We consider the spatially homogeneous but totally anisotropic and non-flat Bianchi type spacetime in the form 

ds 2 = -c 2 dt 2 + A 2 (dx - Kzdyf + B 2 dy 2 + C 2 dz 2 , (2) 

where A, B and C are the scale factors and functions of t only. K is a, constant with dimension length -1 . Choosing 
K = one gets Bianchi type I metric. For studying within the framework of Bianchi type II metric, which we are 
interested here in this paper, it is usual to choose K = 1. We also, for convenience, switch to natural units c = 1 
and87rG = l. 

We assume the overall EMT T? consists of two different components; T^j 3 that shall represent the matter 
fluids (dust, radiation etc.) including Cold Dark Matter (CDM) and T^ e ^ ^ that shall represent the hypothetical 
fluid we consider in the context of DE, i.e., 

T J = 3 + T' dc ) ■? (3) 

We do not expect the fluids to have bulk motion on cosmological scales, hence we assume that the fluids are 
comoving, i.e., fluids' four- velocities are u l = S\ with u l Ui = — 1. Most of the cosmological matter fluids possess 
isotropic pressure and can be represented by a constant EoS parameter. Hence, the EMT that represents matter 
fluid can be written as 

T { f j = diag [-p (m \ p (m \p^\p ( ^} 

= diag [-l,^ m \w (m) ,u;( m V m) , (4) 

where p^ is the energy density, p^ is the isotropic pressure and w^ m ' = constant > is the EoS parameter of 
the matter fluid. However, the nature of the DE is still mysterious and we have neither a priori nor observational 
reasons for confining ourselves to assume DE yields isotropic pressure. Because the fluids are co-moving, one may 
get DE that yields an anisotropic pressure by assuming that the pressure of the DE is proportional to its energy 
density and subsequently the EoS of the DE fluid is direction-dependent. Hence we represent the EMT of the DE 
fluid as follows: 

T«"»' = di ag [- P ^\p^\p^\p^] 

= diag[-l, W i dc ), W ( dc ),u;i dc )]p( dc ) 

= diag [-l,u/ d °) +<5,u/ dc) +7,w (de) ]p (de) ) ( 5 ) 

where p^ is the energy density and p^ e \ Py^ and p{ dc ^ are the directional pressures and wi dc \ and uii dc ^ 

are the directional EoS parameters of the DE fluid respectively along the x, y and z axes. In the third line we 
parametrize the second line by introducing the skewness parameters 5, 7 that determine the deviations of the EoS 
parameters of the DE fluid on the x and y axes from the one on the z axis. Further, because the observations do 
not rule out the possibility of a dynamical and anisotropic DE, we allow all the parameters relevant to DE fluid to 
be time-dependent. 
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In the presence of the EMTs given in (Q} and ([S]), Einstein's field equations (p} corresponding to the metric ^ 
lead to the following set of linearly independent differential equations: 
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B A BC CA B 2 C 2 
where an overdot denotes d/dt. We have a system of five linearly independent equations (jf)j)- (|10[) involving eight 
unknown variables, namely, A, B, C, 

/0 (m) ! ^(do)^ ^(do)^ s 

and 7. Therefore, in order to fully determine the system, 
we need three constraining equations. To do that, we shall make one assumption on the physical ingredients of 
the Universe and two assumptions on the kinematics of the Universe. Before doing so, we need to define some 
kinematical parameters of importance in cosmology. 

The directional Hubble parameters, which express the expansion rates of the Universe respectively along the x, 
y and z axes, can be defined as follows: 



■11,= j, II., ^ and II. ( -. (11) 



The mean Hubble parameter, which expresses the volumetric expansion rate of the Universe, can be given as 

H = ly = l(H x + H y + H z ), (12) 

where V = ABC is the volume scale factor of the Universe. The dimensionless mean DP q, which tells whether the 
Universe exhibits accelerating volumetric expansion or not, is another important kinematical quantity and can be 
defined as 

The Universe exhibits accelerating volumetric expansion if —1 < q < 0, decelerating volumetric expansion if q > 0, 
and exhibits constant-rate volumetric expansion if q = 0. 

The anisotropy parameter of the expansion A is crucial when deciding whether the models approach isotropy 
or not and can be defined as 

2 (T 2 (Hi -H 



A = 3^ = 3 M^fl-l- (14) 



i=l 

where a 2 = ^aija 1 ^ (<xy is the shear tensor) is the shear scalar. A is the measure of the deviation from isotropic 
expansion, such that the Universe expands isotropically if A = 0. Considering the kinematics of the models, we say 
that the model eventually approaches isotropy continuously if i) V — > 00 and ii) A — > as t — > 00 |57| . However, 
considering the physical ingredient of the models, we also expect the energy densities of the each fluids to be positive 
as t — > 00 for a more realistic approach. 

Our first constraint is on the physical ingredients of the Universe. In GR, twice-contracted Bianchi identity tells 
that the overall EMT is conserved, i.e., T 13 ,- = 0; 



P 



^ + 3(1 + w^)p^H + p( d °) + 3(1 + w (dc) )p (dc) H + p {dc ~> (5H X + jH y ) = 0. (15) 
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Regarding the constraint on the physical ingredients of the Universe, it is quite reasonable to assume that the fluids 
are minimally interacting, i.e., the EMTs of the two fluids are conserved separately, as done by the authors of Ref. 
[29] . Hence, for the matter fluid we have 

T (m) ij^ = ^ W + 3(1 + w (m) )p (m) ff = Q) (16) 

which also trivially gives 

p (m) =c {ABC)- {1+w(n,) \ (17) 

where Co is a positive constant of integration. Similarly, from the conservation law for the EMT of the DE fluid we 
have 

T (do) ij^ = ^(do) + 3(1 + w (dc) )j0 (dc )jff + p (dc) ^ 6Hx + = Q) (lg) 

which cannot be solved for p' dc ) trivially as can be done for the matter fluid above. 

For the second constraint (the first kinematical constraint), we consider Bianchi spacetime generalization of 
Berman's special law [5l|6] of variation for the mean Hubble parameter that yields a constant value of mean DP: 

H = £V~*, (19) 

where £ > and n > are constants. This law has been widely used for obtaining exact cosmological models in the 
context of anisotropic DE. This law gives us opportunity to study two kinds of volumetric expansion laws. From 
equations (fT2j) and (fT9|) . we obtain 

V = (net + Cl )i for n ^ 0, (20) 



V = Ae ut 



for n = 0, 



(21) 



where c\ and ci are constants of integration. Thus, the law (|19p provides the power-law (J20J) and the exponential-law 
([2T]) of volumetric expansions of the Universe. One may check that q = n — 1, hence the models for < n < 1, 
i.e., q < can be considered in the context of DE. In this paper, we are particularly interested in the accelerating 
models in the context of DE, hence relatively later times of the Universe rather than the early Universe. Therefore, 
for convenience, we shall only discuss accelerating models and the relatively late time behavior of the Universe in 
detail. Recent observations suggest that today q ~ —0.7 (e.g., see references [58, 59..), hence, assuming that the DP 
changes slowly, particularly the model for n ~ 0.3 would represent the Universe in the vicinity of present time. The 
exponential model (n — 0), on the other hand, would also represent the very late Universe. 

As the third constraint (the second kinematical constraint) to close the system, we assume that the component 
cr 1 1 of the shear tensor <j\ is proportional to the mean Hubble parameter H , following the study of Bali and Jain 
[56j . This constraint leads to the following relation between the scale factors: 



A = (BC) r 



(22) 



where m > is a constant. 

Introducing the three constraint equations, one can now obtain the exact solutions by starting with obtaining 
the scale factor A on the x axis explicitly by using (j2"2"j) and the definition of the mean Hubble parameter (TH?)) with 
([20]) for the power-law solution and with (|2 1 1) for the exponential-law solution. Integrating equation (fTO)) two times, 
we obtain the scale factor B on the y axis in terms of A and V as 



B — c^A exp 



— / —dt\ dt + c 3 f — dt 
V J V J 3 J V 



(23) 



where C3 and C4 are constants of integration. Finally, one can obtain the scale factor C on the z axis by using A 
and B in fl2"2"|) . 

After obtaining all of the scale factors explicitly, one can obtain all other parameters by using the explicit 
functions of the scale factors and the energy density of the matter fluid (fTTj) in equations However, before 

proceeding to obtain the complete exact solutions, here, we would like to discuss two immediate findings. 

First, we would like to emphasize an important point that arises from the second kinematical ansatz, which 
gives (|22[) . and the condition for approaching isotropy continuously. Using the second kinematical ansatz (|22[) we 
obtain 

H x =m{H y + H z ). (24) 
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The condition for approaching isotropy continuously imply that all the directional Hubble parameters must evolve 
to an identical value as t -> oo, i.e., H x = H y = H z at t — oo. One may check from (pH)) that H x = H y = H z can 
be possible only if m = 5. Hence, if m 7^ | one cannot obtain a model that can approach isotropy continuously, 
but a model that approaches isotropy only for a certain period of time and a model that recedes from isotropy 
continuously. However, we also note that m = ^ is the necessary but still insufficient condition for approaching 
isotropy continuously. In other words, choosing m = \ does still not assure approaching to isotropy continuously 
but provides the possibility, i.e., if m — i, one still needs to make further analyses of the anisotropy parameter of 
the expansion whether it can satisfy or cannot satisfy the conditions for approaching isotropy continuously. 

Second, we would like to notice a relation between the BII metric and the anisotropy of the DE fluid. One can 
observe that, from (j6|), (0 and (fT0|) the skewness parameters of the EoS along the x and y axes should be equal; 
i.e., 

7 = 5. (25) 

This shows that Bianchi type II spacetime ([2]) does not allow the pressures of the DE fluid to be different along x 
and y directions as long as the shear stress of the overall EMT is null (i.e., as long as Xf = 0, which gives G 2 = 
(jTUJ) ) . This is a result similar to the one discussed in a recent paper [60J by Saha. Hence, although the metric 
we consider is totally anisotropic, the anisotropy of the DE is constrained to be axially symmetric, as long as the 
overall EMT possesses zero shear stress. 

3 Model for power-law volumetric expansion 

Using the equations (|2H)) . (|2"2"|) and (f2"3")l . we obtain the directional scale factors as follows: 

A = (n^ + ci)^TT) , (26) 

B = c 4 (n£t + ci)'*^™ 1 ) exp 
C = c± x (nit + ci)"(™+i) exp 

where M = (9m — 3 + mn + n)(3m — 3 + mn + n) and n^3. In the special case n = 3, the directional scale factors 
are as follows: 

A = (3£t + Ci )^tt , (29) 



(to + iy 

2£ 2 M 
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C3 



tin - 3) 



(nit + ci)^ 



(27) 



(to + 1)" 
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(28) 



3m( + i 3 (m+l) 

B = c 4 (3tt + ci) 3«<™+i) exp 
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C = c 4 1 (3& + ci) s, i™+ 1 ) exp 



i2m2 (3^ + Cl )^ 
(to + l) 2 



144^2 r 



Using the directional scale factors, we obtain the directional Hubble parameters as follows: 

fl* = -(nit + a)- 1 , 



3m£ , . , 

H v = ^(nit + d)- 1 + 



m+l 
m + l 



(nit + d) »(™+u +1 + C3 ( n £t + ci) 



H z = 



m + l 

3(1 - m)i 



(nit + ay L - 
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m + l 



(nit + a) +1 - c 3 (n& + a) 
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The anisotropy parameter of the expansion is also obtained as follows 
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(30) 
(31) 

(32) 
(33) 
(34) 
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The energy density of the matter fluid is obtained as follows: 



p (m) =co(n£t + Cl y 



(36) 



The energy density, the EoS parameter on the z axis and the skewness of the EoS parameter of the fluid we 
consider in the context of DE, respectively, are obtained as follows: 

M j 9£ 2 m(2 — m) , . . , 15 — 33m — ran — n , „ . e(i-m) 
,(de) = v L(n£ t + Cl y2 + f n £ t + c ^ „ (m+ i) 



(m + l) 2 
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(m + l) 2 
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2c 3 (m + 1) 
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(38) 
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(n£t + Cl y 2 + 2(n£t + Cl ) ^+tj 



(39) 



The density parameter of the matter fluid f2( m ) and the density parameter of DE fluid f2( dc ) may also be given 
as follows: 

2ii-3(l + ra< m > ) 

-(n£t + Cl ) s , (40) 



n (m) _ P_ 



3ff 2 3^ 2 
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+ 2 



(41) 



Adding (|40[) and (|41l) . on defining total density parameter as the sum of the density parameters of the matter 
and DE fluids, we obtain the following: 



Q(T) = fi (m) + fi (de) = 3m(2 m) + (2m l) 2 _ 1 

(m+1) 2 2 



(42) 



One may check that the solution satisfies the EMT conservation equations ([TBI and (118)) . 

The initial time of the Universe is t = —Ci/nl, which can be shifted to t = by choosing c\ — 0. Because this 
study is in the context of DE, we can conveniently consider only relatively later times of the Universe rather than 
early Universe in our discussion about the results. 

We have already discussed that m = \ is the necessary but insufficient condition for approaching isotropy con- 
tinuously and one may, here, check that this model satisfies all the conditions for approaching isotropy continuously 
provided that m = 5 and n < 1 at the same time. In other words, the accelerating models for m = \ approaches 
isotropy continuously and also assures the positivity condition for the energy densities of the fluids in the late Uni- 
verse. We may elaborate the late time behavior of the expansion anisotropy A. One may check that the anisotropy 
of the expansion does not diverge as t — » 00 provided that n < 1 and m < §+^- In our proceeding discussion 
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in this section, we consider only the cases that satisfy these two conditions at the same time, and for sufficiently 
later times of the Universe (for instance, today) we get the following results. Expansion anisotropy converges to a 
constant value for large t values: 

2(2to-1) 2 

(m + \y 

One may easily observe that as would be expected from an accelerating model with m — i, A ps for large t 
values. Otherwise, A converges to a finite non-zero value for large t values. For instance, Campanelli et al. 22 
gives | VA |~ 10~ 5 for the present Universe from WMAP data, which corresponds to m ~ 0.5 ± 0.0000053 in our 
model. 

The behaviors of the energy densities of the fluids also need attention to decide whether our model is realistic 
or not. Obviously for large t values, the density parameter of the matter fluid is 

ps (44) 

independent of whether the Universe exhibits accelerating volumetric expansion or not. One may observe that, on 
the other hand, DE dominates the matter fluid for large t values for accelerating models and its density parameter 
approximates to the following finite value: 

(m + l)2 • 

Hence, the total density parameter also converges to 

(T) 3m(2 - m) 

(m + 1) 2 ■ { ' 

Hence, for sufficiently later times of the Universe, DE energy fluid dominates the matter fluid provided the above 
mentioned conditions, and in particular case for m = 5, we have ps 1. Thus, for m = h and n < 1, our model 
predicts a spatially flat and isotropic Universe at late times. For a demonstration, we depict the anisotropy of the 
expansion in Fig. [I] and density parameters in Fig. [2] In these figures, also in all other figures, we consider m = 5, 
n = 0.27 (q = —0.73) that corresponds to the value obtained from observations for the present Universe (e.g., see 
[SHIES]), w (m) = since the most dominant ingredient after the DE in the present Universe is the pressureless 
matter. On the other hand, for the other constants we have selected the values as i = 2, co = 12, ci = 1, C3 = 1, 
which give the true beginning of the universe t = — However, because we are not interested in the early times 

of the Universe but the later times, for convenience, we plot the graphs from t — to t — 2 that can show us the 
generic behaviors of the parameters for m = | and q = —0.73. From Fig. [TJ one may observe that the anisotropy 
parameter A (|3"5|) has larger values at earlier times of the Universe and it decreases monotonically to zero and 
flattens as t increases (A — > and A — > as t — > 00). One may also check that A diverges at the true beginning of 
the Universe. Moreover, it dominates the physical sources in the sufficiently early times of the Universe. Therefore, 
in view of equation (|42p . the total density parameter (i.e., density parameter of the matter field plus the density 
parameter of the DE) curve goes down as we go back into time (see Fig. [2]). 

We would also like to discuss on the dynamics of the fluid we consider in the context of DE little bit more under 
the conditions n < 1 and m < fe^- The energy density of the matter fluid decreases monotonically as the Universe 
evolves, while the energy density of the DE fluid increases until a certain time of the Universe and eventually 
converges to a finite value depending on the value of m in the later Universe (see Fig. [3]). The directional EoS 
parameters of the DE eventually evolve into the quintessence region. One may observe that both wi de ^ = Wy de ^ and 
w ( dc ) = wi d °^ start in the phantom region and evolve into the quintessence region (see Fig. [¥]). This is consistent 
with the observations, for instance, with the range (—1.45, —0.74) given by Melchiorri et al. [61] . The anisotropy of 
the DE can be measured by w ^ dB) and it asymptotically tends to ^(m+i^-^m? ^m+i) ^ or l ar S e ^ values. In particular, 
(j e) becomes zero at late times for the accelerating models with m = | (see Fig. [5]), which in turn implies that 
the anisotropy of DE goes off during the cosmic evolution. 



^ere and in the figures, for convenience, we have not used dimensions for the cosmic time t, since one can always set the scale and 
the dimension of the cosmic time so as it to correspond to the true cosmic time. 
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Figure 1: The anisotropy of the expansion A versus the cosmic time t. 



1.0 

t 



Figure 2: The density parameters of the dark energy fluid f2( de ) (dotted curve) and matter fluid (dot-dashed 
curve) and the total density parameter fl^ (solid curve) versus the cosmic time t. 



4 Model for exponential volumetric expansion 

Solving the equations ((2T|) , (|22|) and (|23|) . we obtain the directional scale factors as follows: 



Jo?- ( 3ml 

A = c 9 m+1 cxp 1 

2 \m + l 



B = c 4 c^ +1 exp 



3(l-m) 

C = c 2 m+1 exp 



3ml c. 
-t + 



S(m-l) 

m + 1 (™ I ^2 

2 



+ 1 18£ 2 (3m- l)(m- 1) 



(m + 1) 2 c «^ll t _ _£3_ e -3ft 



3^ 



3(1 -to)^ 



6(m-l) 

c 2 m+1 <••• - 1 



m+1 18^ 2 (3to- l)(m- 1) 
The directional Hubble parameters are obtained as follows: 



(m + 1) e m^_ii t + _C3_ e _ m 



3lc\ 



H x 



6(m-l) 
, m + 1 



H,. 



3ml 



m+1 3£(3m-l) 



3ml 
m + 1 



(m + 1) e^(m-i) t C3 



e m+1 



-3ft 



(47) 
(48) 
(49) 

(50) 
(51) 
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Figure 3: The energy densities of the matter fluid pW (dashed curve) and dark energy fluid p^ - 1 (solid curve) 
versus the cosmic time t. 




0.0 0.5 1.0 1.5 2.0 

t 



Figure 4: The directional equation of state (EoS) parameters of the dark energy fluid along the x and y axes 
wi d °^ — Wy (dashed cruve) and along the z axis i(/ dc ) = wi dc ^ (solid cruve) versus the cosmic time t. 



6(m-l) 

_ 31(1 - m) c 2 "' +1 (m + 1) t c z m 

Hz ~~^Tl 3*(3m-l) C T ■ {5) 

The expansion anisotropy parameter is obtained as follows: 

= 2(2m-l) 2 c 2 ' m+l) (5-llm) g^-i) f 4c 3 c 2 "' +1 ' (m + 1) «*(^-') t 
(m + 1) 2 6^ 2 (3m-l) 6 9^ 3 (3m-l) 6 

12(m-l) 

2c 3 (l-2m) 2c 2 m+1 (m+1) 2 i^(^-i) t 2c§ « ft 



^ 2 c|(m + l)" e ' 27£ 4 (3m-l) 2 " ' 3^ 2 



+ ^ffe- M • (53) 



•2 

The energy density of the matter fluid reads as 

p M = c c 2 - 3(1+ "' (m)) e- 3 ^ 1 +«' (n,) ) t , (54) 
The energy density, the EoS parameter on the z axis and the skewness of the EoS parameter of the fluid we 
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Figure 5: Anisotropy of the dark energy fluid S/w^^ versus the cosmic time t 



consider in the context of DE, respectively, are obtained as follows: 

6(m-l) 3(m-3) 

p (dc) = 9^ 2 to(2 - m) + c 2 m+1 (5-llm) c 6€( r+ -i) t _ 2c 3 c 2 m+1 (m + 1) t + 3c 3 (l - 2m) ^_ ut 



(m + 1) 2 4(3m-l) 

12(m-l) 

c 2 m+1 (m + l) 2 i2i(m-i) t c; 



3^(3m- 1) 



C 2 (m + 1) 



9^ 2 (3m- l) 2 



-e m +! 



2 

3,,-6ft ^ -3(lW m) ) -3f(lW m) )t 

fc- OQCq fc- j 



(55) 



w 



(do) 



,(de) 



9^ 2 (m 2 - m + 1) c 2 m+1 (m + 1) ei(m-i) t 2c 3 c 2 m+1 (m + 1) 3< + m - 3 ) t 



(m + 1) 2 4(3m-l) 

12(m-l) 



-e 



3£(3m - 1) 



-e 



3c 3 (l-2m) 3ft c 2 m+1 (m + 1) 2 i2j(m-ii f c 



c|(m + l) 



9^ 2 (3m- l) 2 



e m+1 



c c 2 - 3(1+ ' i,(n,)) e- 3 ^+* ,(nl) ) t 



(do) 



9^ 2 (2m - 1) „ a«m-i) , 
m + 1 ^ 



The density parameters of the matter fluid and the DE fluid are obtained as follows: 



3(l+u> (m) ) 



-3£(l+io (m) )t 



U 2 



(56) 
(57) 

(58) 



q 3m(2 — m) 



(5 — 11m) 6€(m-i) t 2c3C 2 m+1 (m + 1) 3i(m-3) t c 3 (l — 2m) 



(m + 1) 2 12^ 2 (3m-l) 

12(m-l) 

C 2 + 1 (m + l) 2 12<!(m-l) t 



27£ 4 (3m - l) 2 



-e 



-e 



9f 3 (3m- 1) 



e + 



-3ft 



2 c^(m + 1) 



-3(l+tu (m) ) 
C C 2 -3£(lW m) )t 

3^ 2 



The total density parameter is obtained as follows: 

(T) 3m(2- m) + (2m- l) 2 1 
~ (m + 1) 2 2 



A. 



(59) 



(60) 



One may check that the solution satisfies the EMT conservation equations ([TBI and (TT5)) , as expected. 
In this model for n = 0, the DP of the volumetric expansion is q = — 1, which leads to dH/dt = that implies 
the greatest value of Hubble's parameter and the fastest rate of expansion for a Universe expanding forever. This is 
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one of the possible futures of our Universe. Therefore, the derived model can be utilized to describe the dynamics 
of the future of the present Universe. So, in what follows, we emphasize upon the late time behavior of the derived 
model. One may check that only in the models for m < 1, expansion anisotropy A does not diverge as t — > oo but 
converges to a constant, and the energy density of the fluid we consider in the context of DE does not evolve to 
negative values. 

Hence, the detailed picture of sufficiently older Universe is as follows; for m < 1 we obtain 

A* 2(2ro - 1)2 (61) 
(m + 1) 2 ' { ' 

p (m) « 0, (62) 

(do) „ £^m(2-m) (de) _ m 2 - m + 1 (2m-l)(m+l) g _ (2m-l)(m+l) 

~ (m + 1) 2 ' ' ~ m(2-m) ' ' m(2 - m) ' «/ de ) ~ m 2 - m + 1 1 J 

(64) 



(T) _ 3m(2 - m) 
~ (m+1) 2 ' 



and for the particular case m = | , we further obtain 

A w 0, (65) 

p (m) « 0, (66) 
p( do >« 3 * 2 , «,<«*) 5*0, ^y^O (67) 

ft( T ) « 1. (68) 

One may observe that if m = |, as t — > oo, the expansion anisotropy converges to zero, the energy density of the 
matter fluid converges to zero, the fluid we consider in the context of DE mimics the well known positive cosmological 
constant (i.e., DE fluid isotropizes and its energy density takes a positive finite value), the total density parameter 
converges to unity, which means that the Universe becomes spatially isotropic and flat. Thus, the model for m = | 
is well behaved and satisfies all the conditions for approaching isotropy continuously. In short, for m — \ our model 
that exhibits exponential volumetric expansion can evolve to a phase where it is observationally indistinguishable 
from the well known solution of the Einstein field equations in the framework of RW spacetime in the presence of 
positive cosmological constant, which is also known as de Sitter Universe and is the future of the Universe in the 
standard ACDM model. 

5 Concluding Remarks 

We obtain cosmological models that can exhibit rich variety of behaviors in the presence of two minimally interacting 
fluids (a perfect fluid as the matter fluid and a hypothetical anisotropic fluid in the context of DE) within the 
framework of spatially homogenous but totally anisotropic and non-flat Bianchi type II spacetime in GR constrained 
with two kinematical ansatze. As the first kinematical ansatz, we assumed the constant mean DP whose negative 
values give accelerating volumetric expansion. Although it is incapable of representing the whole or a long period of 
history of the Universe at once, it can still be considered as an approximation for representing a particular period 
of the Universe, for instance, the vicinity of the present time and the very late times of the Universe. As the second 
kinematical ansatz, we assumed one of the component of the shear tensor is proportional to the mean Hubble 
parameter. 

We showed that although the Bianchi II spacetime we considered is totally anisotropic, it admits only an axially 
symmetric anisotropic pressure of a comoving fluid as long as shear tension does not exist, independently of the 
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ansatze we considered. We also showed that, for the models that exhibit accelerating volumetric expansion (n < 1), 
there are well-behaved cases where the anisotropy of the expansion converges to finite values (including zero) and 
the fluid we consider in the context of DE yields positive values for its energy density in the late Universe. 

In the particular case with m = i among the well-behaved cases, the anisotropy of the expansion converges 
to zero continuously, the total density parameter converges to unity and the anisotropy of the DE converges to 
zero, which means that our model becomes indistinguishable from a spatially anisotropic and flat Universe in the 
presence of isotropic fluids. Further, in the case for n = and m = 1/2, the fluid we consider in the context of DE 
mimics the cosmological constant and the Universe evolves to the well known de Sitter Universe. 

Yielding anisotropic pressure, the fluid we consider in the context of DE, can produce results that can be produced 
in the presence of isotropic fluid in accordance with the ACDM cosmology. However, our model gives additional 
opportunities by being able to allow kinematics that cannot be produced in the presence of fluids that yield only 
isotropic pressure. There are also well behaved cases where the anisotropy of the expansion and the anisotropy of 
the fluid do not cease but converge to non-zero finite values in the late Universe. Hence, the dominance of such 
anisotropic DE candidates can distort the isotropy of the Universe, which is believed to be gained in the inflationary 
era in the early Universe, with required limits and characteristics. For instance, recently, Campanelli et al. [2"2"] 
investigated the possibility that an anisotropic DE component could generate an ellipsoidal Universe with the correct 
characteristics to explain the low quadrupole in the CMB fluctuations and proposed that an ellipsoidal Universe 
allowed better matches with the large-scale CMB anisotropy data. They give | V~A |~ 10~ 5 for the anisotropy of 
the expansion of the present Universe from WMAP data, which corresponds tom~ 0.5 ± 0.0000053 in our model. 

Our models where the Universe accelerates with an expansion anisotropy that can evolve to a finite value 
(including zero) at late times would be of interest to many readers since it can generate isotropic expansion as 
would be expected in standard cosmology as well as the non-divergent anisotropic expansion in accordance with 
the discussions on whether isotropic expansion is good enough to explain cosmological observations or not. 
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